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ABSTRACT 

We analyze a large collection of RXTE archive data from April 1997 to August 2003 of the bright X- 
ray source Scorpius X-1 in order to study the broadband spectral evolution of the source for different 
values of the inferred mass accretion rate by studying energy spectra from selected regions in the 
Z-track of its Color-Color Diagram. A two-component model, consisting of a soft thermal component 
interpreted as thermal emission from an accretion disk and a thermal Comptonization component, is 
unable to fit the whole 3-200 keV energy spectrum at low accretion rates. Strong residuals in the 
highest energy band of the spectrum require the addition of a third component that can be fitted with 
a power-law component, that could represent a second thermal Comptonization from a much hotter 
plasma, or a hybrid thermal/non-thermal Comptonization, where the electrons in the Comptonizing 
cloud have a Maxwellian distribution of velocities with a power-law hard tail. The presence of this 
hard emission in Scorpius X-1 has been previously reported, however, without a clear relation with 
the accretion rate. We show, for the first time, that there exists a common trend in the spectral 
evolution of the source, where the spectral parameters change in correlation with the position of the 
source in the CD. In particular, using a hybrid thermal/non-thermal Comptonization model (EQPAIR 
code), we show that the power supplied to the non-thermal distribution can be as high as half of the 
total hard power injected in heating the electron distribution. We also found that a small sample of 
spectra, when the source resides at the top of the FB can also show intense hard X-ray emission. We 
discuss the physical implications derived from the results of our analysis, with a particular emphasis 
on the hardest part of the X-ray emission and its possible origins. 

Subject headings: accretion discs - stars: individual: Scorpius X-1 — stars: neutron stars — X-ray: 
stars — X-ray: general — X-ray: binaries 
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1. INTRODUCTION 

1.1. Scorpius X-1 

Scorpius X-1 is the brightest persistent X-ray source 
in the sky and the f irst identified X-ray extra-solar 
source (jCiacconi et al.l 1 1962). The X-ray source is an 
old, low magnetized neutron star (NS), accreting mat- 
ter transferred through Roche-lobe overflow from a low- 
mass companio n (recently identified as a n M class star 
of - 0.4 M©, ISteeehs fc CasaresI [200l) . Scorpius X- 

1 is a prototype of the class of the Low-Mass X-ray 
Binar ies (LMXBs), and as suming a distance of 2.8 ± 0.3 
kpc (Bradshaw et al.lll999| ). the source emits close to the 
theoretical Eddington limit for a 1.4 Mq NS {LEdd ^ 

2 X 10^^ erg s'^). 

Based on the timing behavior of LMXBs in correla- 
tion with the position of a given source in the X- 
ray co lor-color diagram (CD), iHasinger fc van der Klid 
(|1989D grouped these sources into two categories: the Z 
sources and the Atoll sources. The former are brighter, 
radiating close to Lsdd, the latter are less bright, emit- 
ting at 0.01-0.1 LEdd- Z sources exhibit the classical 
three-shaped branches describing a Z pattern in the CD: 
the Horizontal Branch (HB) at the top of the Z track, 
followed by the Normal Branch (NB) and the Flaring 
Branch (FB) at the bot tom of the pattern . There are 
strong indications (e.g. IVrtilek et al.l Il991h suggesting 

^ Dipartimento di Scienze Fisiche ed Astronomiche, Univer- 
sita di Palermo, Archirafi 36, 90123 Palermo, Italy; email: 
dai@fisica.unipa.it 

^ Nicolaus Copernicus Astronomical Center, Bartycka 18, 00-716 
Warsaw, Poland 



that what drives the changing in the spectral and tempo- 
ral properties of LMXBs is the instantaneous accretion 
rate (M), which, for Z sources, is believed to increase 
monotonically from the HB to the FB. Similarly, Atoll 
sources, display two different spectral/timing states: the 
soft and luminous banana state, associated to higher 
accretion rates, and the hard and less luminous island 
state, associated to lower M. This early and straight- 
forward classification has now moved to a more complex 
picture, as the patterns that Atoll and Z sources describe 
in their CDs, when displayed on long timesc ales, appear 
to be more and more similar t o each other ()Muno et al.l 
[200l iGierlihski fc Donel[200l ). although some observa- 
tional facts still hold: Z sources move on the CDs on 
shorter timescales, have higher luminosities and present 
a timing phenomenology different from that of the Atoll 
sources. 

Scorpius X-1, as all the sources of its class, shows several 
quasi-periodic oscillations (QPOs) along all the branches 
of its CD: the horizontal branch oscillations (HBOs), 
the normal branch oscillations (NBOs) and the flaring 
branch oscillations (FBOs). In Scorpius X-1 NBOs, with 
peak frequencies in the range 4.5-7 Hz, and FBOs, in the 
frequency range 6-25 Hz, seem to be physically related 
to each other since the NBO peak frequency smoothly 
joins the FBO peak frequency as the so urce moves from 
the NB to the FB (jCasella et al.l 120061 ). Van der Klis 
et al. (1996) reported the first observation of HBOs in 
Scorpius X-1 at 45 Hz, with an inferred harmonic near 
90 Hz; the power spectrum can also show a pair of QPOs, 
whose frequencies are in the range 800-1100 Hz, denoted 
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kHzQPOs, that shift simultaneously in frequency, with 
an almost constant or weakly frequenc y-dependent peak 
separation (see e.g. lZhang et al.|[2006[ ). The physical in- 
terpretation of these timing features is no t unique and 
the r elated scientific debate is still open (see Ivan der KlisI 
l2004l for a review). 

Spectral studies of Scorpius X-1 have been so far not 
so extensive and detailed as timing studies. This is in 
part due to the strong brightness of the source that ac- 
tually prevents its observation with the most sophisti- 
cated and high-resolution X-ray satellites like BeppoSAX 
and ASCA in the recent past and Chandra and XMM- 
Newton at present. Three articles have investigated so 
far the broad-band spectral behavior of the source, us- 
ing data from the Ro ssi X-ray Timing Explorer (RXTE) . 
iBarnard et al.l ()2003[ ) used data from both the Propor- 
tional Counter Array (PCA) and the High Energy X-ray 
Timing Experiment (HEXTE), showing that the spec- 
trum, in the energy range 2.5-50 keV, can be fitted by a 
simple two-component model, a black-body soft compo- 
nent and a cutoff power-law, plus a broad Gaussian line, 
and interpreted these resu lts in the framework of the so 
called Birmingham model (jChurch &: Balucihska-ChurchI 
1200 If) . In this interpretation the spectral emission con- 
sists of black-body emission from the NS plus Comp- 
tonized emission from an extended accretion disk corona. 
iBradshaw et all ()2003D used only PCA data (in the 2.5- 
18.2 keV range), adopting a model consisting of a black 
body emission plus a bulk motion Comptonization com- 
ponent and a reflection component in the form of a broad 
Gaussian line. D'Amico ct al. (2001) studied HEXTE 
data above 20 keV, in order to test the presence of high 
energy X-ray emission in the spectra of the source. Data 
were fitted, using a bremsstrahlung component, to mimic 
the effects of the thermal component, and a power-law 
component, whenever the statistics and clear residuals in 
the fit required it. In this way no apparent correlation 
was found between the presence of the hard tail and po- 
sition of the source in the CD, contrarily to what had 
been reported for other Z sources. 

More recently, Di Salvo ct al. (2006), using INTEGRAL 
data in the 20-200 keV energy band, detected a hard 
X-ray tail, with photon index values between 2 and 3.4, 
whose intensity decreased as the source moved towards 
higher accretion rates. Data did not show evidence of a 
high-energy cutoff up to ~ 200 keV, suggesting a non- 
thermal origin for the hard tail. 

Scorpius X-1 is also an inter esting source of radio 
emission. iFomalont et~all (|2001f ). through an extensive 
VLBA monitoring campaign, have shown that the radio 
emission is composed of a point-like radio core emission 
at the position of the X-ray source and of two opposite 
radio lobes, moving through the ISM with relativistic 
speeds v/c = 0.45 ± 0.03, and with an angle of 46°±6° 
with respect to our line of sight. It is also evident a 
connection between the phenomenon of radio flare at the 
core of the source and the following flaring of the lobes, 
so that it is argued that the energy production of the 
radio emission is confined near the NS and only after- 
ward transported to the lobes via the working surface of 
a jet (see e.g. the magnetohydrodynamical simulations 
of kato et ahllMil . 

These results stimulated the search for the possible con- 
nection between inflow and outflow mechanisms work- 



ing in the violent regimes of accretion onto a NS or 
a black hole (BH). Energetic mass outflows, collimated 
sometimes in a typical jet pattern, can be produced not 
only in extragalactic X-ray sources, such as quasars and 
active galactic nuclei (A GN), but al so in galactic X- 
ray binary systems (see iFended [20021 for a review). It 
has been observed that galactic BH binary systems are 
able to produce strong jet emission (sources for which 
a jet has been already spatially resolved are usually de- 
noted microquasars) , with radio-loud states associated 
to hard/low X-ray states, alternated to periods of radio- 
quenching during soft/high X-ray states. For NS sys- 
tems radio-loud episodes are generally interpreted as a 
jet signature. In this context, all the Z sources seem to 
owe, under certain conditions, a radio-jet nature. De- 
tections of radio-loud states are usually associated with 
spectral states of low M , i.e. on the HB of the Z-track 
(jFender fc Hendrvll2000( l. while the few radio detections 
in atoll systems fo und the sources resid ing in the island 
state of their CDs (jMighari et al.ll2003[ ). 

1.2. Spectral properties of Z sources class 

We rapidly review in this section the spectral proper- 
ties of the Z sources. The class of the Z sources display 
a homogeneous pattern of spectral properties; their 
energy spectra can be usually decomposed as follows: 

• a soft thermal component, with characteristic tem- 
peratures in the 0.5-1 keV range, interpreted as 
thermal emission from the NS surface or an opti- 
cally thick, geometrically thin accretion disk; 

• a thermal Comptonized component, where the elec- 
tron cloud distribution has a temperature in the 
2.5-3 keV range, the soft seed photon temperature 
is around 1 keV and the optical depth can have 
rather large values (r > 3); this relatively cold ther- 
mal Comptonization is thought to take place close 
to the NS surface; the soft seed photon temperature 
generally exceeds the highest temperature reached 
in the accretion disk, so that the source of the 
soft emission is probably confined in the boundary 
layer between the inner edge of the accretion disk 
and the NS surface. The high optical depth val- 
ues associated to the cloud saturates the seed pho- 
ton spectrum close to the electron temperature, so 
that a quasi Wien spectrum results. A theoretical 
interpretation of this spectral decomposit ion has 
been given by llnogamov k SunvaevI (jl999| ). where 
the boundary layer emission is expected to be ra- 
diation pressure supported and its emission is lo- 
cally at the Eddington rate. Change in luminos- 
ity is attributed to a change in the emitting area 
rather than in its emissivity. Phase resolved spec- 
troscopy of bright LMXBs (Revnivtsev & Gilfano^ 
l2005l ) further supports this interpretation, estab- 
lishing in all the analyzed Z sources a common cut- 
off, due to the boundary layer thermal emission, at 
~ 2.4 keV. 

• a reflection component, often simply modelled with 
a broad Gaussian line, in the 6.4-6.9 keV range. 
In all the Z sources this line has been always ob- 
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served; spectra from BeppoSAX or RXTE observa- 
tions were not able to constraint the shape of the 
line, while the relative width was in the range 0.1-1 
keV. High resolution spectra in this energy range 
have been ob tained so far only with the BBRXT 
(ISmale et all 11993.) an d the XMM-Newton satel- 
lites ()Costantini et all [2005?) for Cy gnus X-2; in 
both observations the line was associated to highly 
ionized iron (centroid energy ~ 6.7 keV), while the 
line a ppears to be intr insically broad (FWHM ^ 1 
keV). iBrandt fc Matt! (1994) pointed out that the 
determination of the line profile can be a primary 
diagnostic tool, but, at the same time, kinematic 
and relativistic effects can greatly distort the line 
from the simple Gaussian profile. Moreover, if re- 
flection is caused by hard X-rays reflected by a cold 
accretion disk, the contribution of the Compton re- 
flected continuum should be taken into account, 
when the iron line emission is a considerable con- 
tribution to the total energy flux; 

• a power-law hard tail, which can contribute up to 
few percent to the total energy flux, whose strenght 
usually varies in correlation with the position of the 
source in the CD, namely being the strongest on the 
HB, gradually decreasing as the source moves to 
higher accretion rates and totally fading in the FB 
(e.g. GX 17+2, Dl_Salvo_et al. 2000). The photon 
index of the power-law was generally found in the 
range 1.5-3 with no evident high energy cutoff up 
to energies of about 200 keV. 

In this article we report a complete investigation of the 
spectral properties of Scorpius X-1 through an extensive 
analysis of RXTE archive data, indicating a clear con- 
nection between position of the source on the CD and 
spectral behavior. 

2. DATA REDUCTION AND ANALYSIS 

The scientific payload of RXTE consists of three in- 
struments, the PCA, the HEXTE, and the All Sky Mon- 
itor (ASM). The PCA consists of five co-aligned Xenon 
proportional counter units (PCUs) with a total effec- 
tive area of about 6500 cm^. The instrument is sen- 
sitive in the energy range from 2 keV up to 60 kcV 
(jJahoda fc PCA Team|[l996r ) , although the response ma- 
trix is best calibrated in the energy range 3-22 keV. 
Data can be processed using several different configura- 
tion modes; for our analysis we exclusively use the Stan- 
dard2 mode, with 16 s time resolution and 128 energy 
channels in the 2-60 keV energy range. The HEXTE 
consists of two cluster of four Nal/CsI-phoswich scintil- 
latio n counters that are sen sitive from 15 keV up to 220 
keV (jRothschild et al.lll998f) . We use the Standard Mode, 
with 64 energy channels, for the reduction and analysis 
of the HEXTE data. Background subtraction is done by 
using the source-background rocking of the collimators. 
We use HEXTE response matrices of 1999 August. 
We collected a large amount of RXTE archive data, dis- 
carding only minor shorter observations, from April 1997 
up to August 2003. We present in Table[T]the datasets we 
used for our analysis, indicating the associated proposal 
number, the starting and ending times of each observa- 
tion in Terrestrial Time, and the corresponding exposure 
times. 



Data have been processed using the standard selection 
criteria, discarding data taken during Earth occultations 
and passages through the South Atlantic Anomaly. We 
only used data from PCU2 for the PCA and data from 
Cluster A for the HEXTE instrument. We constructed 
color-color diagrams (CDs) of the source by extracting 
energy-dependent lightcurves using PCA energy chan- 
nels 5-10, 11-16,17-22 and 23-44, with a 64 s bintime. 
These channel ranges correspond to the energy ranges 
1.94-4.05 keV, 4.05-6.18 keV, 6.18-8.32 keV, 8.32-16.26 
keV respectively. We define the Soft Color (hereafter 
SC) as the ratio of count rates in the 4.05-6.18 keV and 
1.94-4.05 energy bands, while the Hard Color (hereafter 
HC) as the ratio of count rates in the 8.32-16.26 keV 
and 6.18-8.32 keV energy ranges. However, the channel 
to energy conversion depend on the period of activity of 
the satellite (there are 5 different instrumental Epochs), 
most of our data belong to Epoch 3 for which the given 
energy bands are appropriate; for the Epoch 5 datasets 
the energy ranges are slightly shifted, so that this re- 
sults in a consequent shift on the HC and on the SC 
axes. Because we are not interested in the secular shifts 
of the track in the CD, nor we want to accumulate spec- 
tra taken at different Epochs, but to perform a statistical 
study of the broad-band spectral behavior of the source 
for different CD positions in different datasets, we did 
not perform the energy dependent corrections to the col- 
ors. 

We selected regions in the CDs for each observation 
dataset, in order to cover a homogeneous part of the Z- 
track and to have at the same time a suitable statistics. 
From these selections we derived the good time intervals 
(GTI) which we used for extracting the corresponding 
spectra for the PCA and HEXTE data. 
Given the high luminosity of the source all the PCA spec- 
tra have been dead-time corrected. Some observations 
were performed with the use of an offset pointing; in 
these cases we extracted the responses matrixes of both 
the HEXTE and the PCA instruments, following the in- 
dications given in the on-line Data Analysis documenta- 
tion pages. We processed and analyzed data using ver- 
sion 6.0 of the FTOOLS package suite and version 11.3.1 
of Xspec. 

3. THE DATA 

Hereafter we shall refer to a collection of close in time 
pointed observations used to extract one CD by using 
its proposal number. From datasets 20035 (and 30035) 
we have obtained two different CDs; we will label them, 
respectively, 20035A and 20035B (30035A and 30035B, 
see the Spectrum column in Tabic [2]). Datasets 20035A, 
20035B, 30036, 30035A, 30035B and 40020 belong to 
Epoch 3, while datasets 70014, 70015 and 80021 belong 
to Epoch 5. Datasets 70014 and 70015 have been merged 
in one CD, because these observations were close in time 
and the source pointing was the same. 
The eight CDs that we have extracted present the source 
mostly on the NB/FB, as these are the states where the 
source spends most of the time. In Scorpius X-1 the 
HB is not horizontal but rather a vertical continuation 
of the NB (a characteristic shared also by GX 17+2 and 
GX 349+2, the so-called Sco-like sources). For visual 
clarity we will refer to the Scorpius X-1 CD pattern, 
mainly making the distinction between the left (or FB) 
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TABLE 1 

LOG OF THE RXTE OBSERVATIONS USED IN THE ANALYSIS 



Proposal 


Starting time (TT) 


Ending time (TT) 


Effective Time 


Number (CD) 






ks 


20035A 


1997-04-18 15:20:32 


1997-04-24 19:02:40 


75.6 


20035B 


1998-01-03 21:17:52 


1998-01-08 13:19:44 


60.0 


30036 


1998-01-07 13:29:52 


1998-01-08 22:16:48 


30.4 


30035A 


1998-05-30 00:28:00 


1998-06-02 07:01:36 


47.8 


30035B 


1998-07-02 06:26:40 


1998-07-05 10:32:00 


32.2 


40020 


1999-01-06 05:25:20 


1999-01-16 23:34:40 


122.1 


70014 


2002-03-15 01:02:56 


2002-05-16 23:20:32 


139.0 


70015 


2002-05-19 21:55:44 


2002-05-19 23:54:40 


2.24 


80021 


2003-07-30 03:00:48 


2003-08-13 16:43:44 


146.9 



Note. — Observation log of the RXTE datasets used for our analysis. Some datasets, although having the same proposal number 
contain pointed observation separated by long time intervals. In these cases we collected all the pointing observations close in time in just 
one CD, and distinguish the relative CDs by a letter, following the proposal number. Starting and ending times are expressed in Terrestrial 
Time. Effective exposures times are filtered by our GTI screening criteria. 



track and the right (NB and HB) track, and between the 
top (harder) parts and the bottom (softer) parts of the 
tracks. 

The count rates associated with the extracted spectra 
show a significant trend, generally correlated with the 
HC; HEXTE and PCA count rates are the highest at the 
top of the FB and monotonically decrease towards the 
FB/NB apex; following the direction bottom NB — » HB, 
the trend is inverted. This behavior is mostly stressed in 
the HEXTE data (compare the countrate in spectrum 07 
in the top FB with spectrum 10 in the transition region 
FB/NB, for which the rate is reduced by two thirds in 
the HEXTE, but only by one third in the PCA). 

In four datasets (20035A, 20035B, 40020 and 30036) 
the satellite pointed the source on axis; in three datasets 
(30035A, 30035B, 80021) an offset pointing has been 
used. This influenced the high energy statistics, where in 
the same position on the CD, an offset pointing produced 
an HEXTE count rate lowered by almost one order of 
magnitude. Spectra extracted from datasets 70014 and 
70015 have also low statistics in the HEXTE band, as 
the integration time is shorter compared with the other 
spectra. 

We generally used a 3.0-22 keV energy band for the PCA 
spectra and a 20-200 keV energy band for the HEXTE 
spectra. HEXTE channels from 28 to 63 have been re- 
binned, grouping four channels into one. In some of the 
spectra that we extracted we found, however, model- 
independent mismatches in the overlapping energy re- 
gion between the HEXTE and the PCA; because these 
mismatches are of instrumental origin we occasionally 
restrict these energy bands until these features did no 
longer impact on our fit results. In Table [2] we report the 
exact energy bands used for the PCA and the HEXTE 
datasets for the selected spectra. 

The choose to assign an appropriated systematic error 
for the PCA data is an essential step in the analysis of 
the spectra. This operation is not straightforward, as 
the commonly used standard candle, the Crab, provides 
in this case only a partial calibration tool. The Crab 
countrate for every energy channel below ^ 20 keV is 
about one order of magnitude less than the Scorpius X- 
1 countrate, thus making the extrapolation of the Crab 
based systematic error to the Scorpius X-1 spectra risk- 
ful. The mostly adopted view to assign a 1% system- 



atic error for all the energy channels has a strong impact 
on the softer part of the spectrum, where there are the 
mostly weighted channels; we found that in this case the 
fit is strongly driven by residuals in the 3-6 keV range, 
leading to a modelization of the spectrum that does not 
seem physically plausible. The overall effect is a shift of 
any soft component outside the RXTE energy band, and 
unrealistic high values for the extrapolated flux. 
The only way to avoid the strong driving effect of the 
first channels is to arbitrarily raise the systematics in 
this range; we found that assigning a 3% systematic er- 
ror for channels between 3 and 3.9 keV (namely the first 
two PCA channels available in our spectra), while leav- 
ing a 1 % systematic error for all the other channels, al- 
lows to derive a physically plausible scenario for the spec- 
tral evolution of the source, does not result in extreme 
lower/higher values of reduced for every fit performed, 
allows us to exploit all the available energy band and to 
generally constrain all the spectral components of the 
adopted models. We occasionally still found in some fits 
systematical residuals in the soft part of the spectrum, 
that are more stressed in the softer spectra, whose origin 
can be presumedly instrumental and related to the xenon 
L-edge near 5 keV. We choose to ignore these effects, in 
order to not overparameterize the fit, as they do not in- 
fluence the determination of the spectral parameters and 
do not impact the values of the fits. 
No systematic error has been associated to the HEXTE 
data. A normalization constant is left free to vary be- 
tween the PCU2 and HEXTE spectra to take into ac- 
count residual flux calibration uncertainties. 

4. SPECTRAL MODELS 

As pointed out in Section 1, spectra of LMXBs are 
usually described as the sum of at least two spectral 
components. Because the PCA energy band starts from 
3 keV, it is not possible to constrain the effect of the 
photoelectric i nterstellar abso r ption on the source flux. 
Following Christian fc Swanld (|1997t l. we fixed its value 
to 3 X 10^^ cm~^, for each fit performed. To fit the CDs 
selected spectra, we tried at first a series of models given 
by the combination of a soft thermal component, such 
as a bl ackbody or multicolo r disk blackbody (DISKBB in 
Xspec. lMitsuda et al.lll98^ . and a thermal Comptonized 
component. For the latter we, alternatively, tried the 
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Fig. 1. — CDs extracted from RXTE datasets used in our spectral analysis. Different colors indicate different selected regions 
from which the CD-resolved spectra were produced. We label the spectra through a progressive numbering as shown in the 
CDs. 

of the power dissipated along the disk, if t he disk is very 
close to the NS (see also lDone et al.|[2002l for a more ex- 
tensive discussion). 

The DISKBB has two free parameters, the temperature 
expressed in keV at the inner disk radius {HTdb) and a 
normalization factor that depends on the inner radius of 
the accretion disk, the distance from the source and the 
inclination angle between the line of sight and the normal 
to the disk; the COMPTT has four free parameters: the 
soft seed photon temperature, kTo, the electron temper- 
ature of the Comptonizing cloud, kTg, the optical depth, 
r, and a normalization constant. We assumed that the 
Comptonizing geometry is spherical. We will refer to this 
two-component model as the DBBTT model. 
As it is evident from Table ^ the DBBTT model is, 
however, gradually unable to satisfactorily fit the whole 
energy spectrum as the source resides in zones of low 
inferred accretion rate. The residuals, with respect to 
this model, clearly indicated an excess of flux at ener- 
gies greater than 30 keV. There is also a second group of 
spectra (namely spectra 7,8,12,f3 and 26) that gave for 
the same reason an unsatisfactory high value. These 
spectra lie at the top of the FB in the datasets that have 
the higher HEXTE statistics. We will refer to this par- 
ticular set as the top-FB spectra. 

Because the mechanisms producing this hard excess in 
the Z sources are not clear yet, in the following we will 
adopt two working hypothesises to fit the broadband con- 
tinuum of the source: a) as in the case of the quasi satu- 
rated boundary layer Comptonization, this component is 
also the result of a thermal Comptonization process, but 
from a much hotter plasma; b) a broad band spectrum 
is the result of a hybrid electron velocity distribution, 
i.e. a thermal electron distribution, which produces the 
optically thick Comptonization, with a power-law non- 
thermal tail which is responsible for the high-energy ex- 
cess. Because, for case a), the hot plasma can only re- 
side in regions close to the central compact object, where 
magnetic and gravitational effects are the strongest, we 
assume that the boundary layer radiation field is the 
primary source of soft-seed photons also for this second 



COMPPS fPouta nen et al.lll996D . COMPTT ()Titarchuklll994D 
and THCO MP (not included in t he standard XSPEC pack- 
age, see IZdziarski et al.l[l996l ) Comptonization models. 
For any model adopted, large residuals in the 6-10 keV 
energy range indicated the presence of iron features, that 
we simply modelled with a broad Gaussian line in the 6- 
7 keV energy range. We noted a general tendency for 
any adopted model to have a very broad Gaussian line 
(with a > 0.8 keV) and, generally, with centroids at en- 
ergies < 6.4 keV. In order to avoid unrealistically large 
Gaussian lines, absorbing the underneath continuum, we 
constrained the line to have centroid energies only in the 
6.4-7.0 keV energy range, with a line width less than 0.8 
keV, while the normalization of the line was a completely 
free to vary parameter. 

We found an adequate description of the data, for en- 
ergies below 30 keV, using a two-component model 
given by the combination of the DISKBB and the 
COMPTTcom ponents, a modelizati on that was firstly in- 
troduced bv lMitsuda et al.1 (|1984f) . and that is known as 
the Eastern model. Other Comptonization models can 
also statistically equally well represent the spectral sam- 
ple, but we preferred to illustrate our results using this 
Compton model as it is one of the mostly used in the 
LMXBs spectral analysis and allows an easy way to com- 
pare our results with results obtained for other sources. 
Using t he other often use d modelization, the Western 
model (jWhite et all Il986f) that predicts a ~ 1-2 keV 
blackbody thermal emission and a Comptonization com- 
ponent with disk s oft seed-photon in th e 0.4-0.8 keV en- 
ergy range (see e.g lBarnard et al.ll2003D . we found higher 
values for each spectrum and an extrapolated super- 
Eddington luminosity along each part of the Z-track, that 
would be dramatic for some FB spectra which would 
overcome by more that an order of magnitude the Ed- 
dington limit, thus making this assumption not realistic. 
Moreover, the power dissipated in the Comptonizing disk 
corona would be constantly much greater than the power 
dissipated as thermal emission on the NS-surface and this 
breaks our expectations that the power dissipated near 
the compact object should be at least of the same order 
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thermal Comptonization. We note that this assumption 
is also needed in order to avoid an overestimation of the 
contribution of steep hard tail at low energies. The use of 
an unbroken power-law at steep photon indeces can give 
a major contribution in the lowest energy band of the 
spectrum so that the fit is driven by low-energy residuals, 
possibly resulting in an erroneous continuum modeliza- 
tion. Thereafter, a low-energy cutoff at the seed-photon 
temperature must be taken into account as the expected 
photon field spectrum is within the instrumental range. 
We added, thereafter, to our basic two-component model, 
a power-law in order to mimic this second Comptoniza- 
tion, multiplied by a low-energy cut-off. We set the 
value of the low-energy cut-off anchored to the kTg value 
of the COMPTT component, while we choose the power- 
law to have a pegged normalization expressed as flux 
of this component in the 20.0-200.0 keV energy range 
(PEGPWRLW, in XSPEC). We checked for each spectrum 
the improvement by adding a high energy cutoff, not- 
ing that we could obtain only lower limits on the cutoff 
energy, generally in the 30-100 keV range without any 
significant improvement in the x^. Thereafter we did not 
use a high energy cutoff. We will refer at this model sim- 
ply as DBBTTPEG or two Comptonizations model. 
To model the possibility given in case b) we also fitted 
the spectra showing a strong hard excess (namely spec- 
tra for which the DBBTT model gave a null hypotesis 
probability less than 0.05 (it generally corresponded to 
a reduced x^ value of 1.4) with a model given by the 
sum of a thermal disk component plus the recently devel- 
oped thermal/nonth ermal hybrid Comp tonization model 
named EQPAIR (see ICoppilllQQOl 120001 for a fuU descrip- 
tion of the model). It embodies Compton scattering, 
pair production and annihilation, pe^ and e^e^ ther- 
mal and non-thermal bremsstrahlung, and energy ex- 
change between thermal and non-thermal parts of the 
distribution via Coulomb scattering. Selected elec- 
trons are accelerated to suprathermal energies and the 
thermal part of the distribution can be additionally 
heated. This model assumes a spherical plasma cloud 
with isotropic and homogeneous distribution of photons 
and , and soft seed photons produced uniformly within 
the plasma, with a thermal temperature of kTsB- The 
properties of the plasma depend on its compactness, 
/ = Cot I {'R-TTieC^) where £ is the total power in the 
source, TZ is the radius of the sphere and ctt is the Thom- 
son cross section. The compactness is divided in a hard 
compactness, Ih, which corresponds to the power sup- 
plied to the electrons, and a soft compactness, 4, corre- 
sponding to the power supplied in the form of soft pho- 
tons. The compactnesses corresponding to the electron 
acceleration and to the additional heating of the ther- 
mal part of the distribution are denoted as Inth and 
Ith, respectively, and 4 = Inth + kh- The non-thermal 
energy distribution of the electrons in the plasma is as- 
sumed to be a power law replacing the Maxwellian ex- 
ponential tail from a gamma Lorentz factor = 1.3 
up to a ^max — 1000 with an energy index Ginj- The 
total Thomson optical depth and the electron temper- 
ature of the plasma are computed self-consistently from 
the assumed optical depth Tp of the background electron- 
proton plasma. 

We choose to model the thermal disk emission using as 
in the hot Comptonization model the DISKBB compo- 



nent for spectra in HB/NB , while we adopted th e post- 
Newtonian model DISKPN (|GierHnski et al.lll999f ) in the 
case of the topFB spectra, as, in this accretion state, we 
expect the disk to reach its minimum inner disk radius, 
close to the NS; the DISKPN model takes into account 
the effective gravitational potential in the neighborhood 
of the compact object by computing the emergent disk 
spectrum using a more appropriate post-Newtonian po- 
tential. We assumed that the disk reaches the last sta- 
ble orbit at 6 Rg and fixed this parameter for the fits 
performed, while we left free to vary the maximum disk 
temperature {T^max) and a normalization factor (K = 
]VP cos{i) / D"^ (3'^ , where M is the central mass expressed 
in solar units, i is the inclination of the disk, D is the 
distance to the source in kpc, (3 the color/effective tem- 
perature). We have left free to vary for the EQPAIR com- 
ponent the following parameters: kTsB, Ih/h, Inth/lh, 
Tp and the normalization factor. The Ginj was not very 
sensible to variations in the range between 1.5 and 3 for 
spectra in the HB/NB, while for the topFB spectra had 
to be left as a free to vary parameter, as this gave a sig- 
nificant decrease in the when left as a free to vary 
parameter. The value of the soft compactness 4 was un- 
constrained by the fit (it is mainly driven by the pair 
production rate, which manifests itself with the annihi- 
lation line at energies ~ 500 keV, obviously far beyond 
our energy limits), and we kept i t frozen at the value of 

10 (see also iGierlihski et al.lll999t ). All the other param- 
eters of the model were frozen at the default values. We 
will refer to this modelization as the DBBEQPAIR model, 
or simply, the hybrid Comptonization model. 

5. RESULTS AND DISCUSSION 

We have examined 43 energy spectra of Scorpius 
X-1 in the 3.0-200 keV energy band using RXTE 
pointed observations from April 1997 to August 2003; 
the spectra have been extracted from selected regions 
chosen in the X-ray CD. We produced a CD for each 
close in time dataset in order to avoid shifts of instru- 
mental origin, and repeated our analysis for each CD 
obtained in this way, thus having at hand a robust 
representation of the spectral evolution of the source in 
all its accretion states. 

To fit the spectra we firstly adopted a two-component 
model given by the sum of a soft thermal disk emission 
and a hard Comptonizcd emission. We checked for the 
presence of harder X-ray emission by adding a high- 
energy power-law to our basic model, while for a subset 
of spectra we adopted a hybrid thermal/non-thermal 
Comptonization as described in the previous section. 
We find that the basic two-component model DBBTT can 
fit the 3-200 keV energy band spectrum (fits for which 
we obtained a xled — l-^) in 25 out of a total 43 spectra; 
these spectra mostly belong to the FB and to the NB. 
Hard tails, in this group of spectra, are not significantly 
detected in 14 spectra (these are spectra 1, 9, 14, 15, 17, 
18, 22, 23, 24, 29, 34, 35, 39, 40, see Table [l]) when the 
source lies in the FB/NB apex of its CD. For the other 

11 spectra, although the addition of a new component 
is not strictly required by the fit, we found a general 
improvement in the value. 

On the contrary, we find that this model fails to describe 
the spectra in 18 of the 43 selected spectra (fits for 
which we obtained a xled — l-^); t^ie fits generally 
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TABLE 2 

COLOR-DIAGRAM RESOLVED SPECTRA AND FITS RESULTS 



Q 

opectrum 


Exposure 


Countrates (cts/s) 


Energy 


range (^KeVj 


HRRTT 




nRRpriD A TR 
UDDlLmrii. In, 




s 


PCA 


HEXTE 


PCA 


HEX IE 


Xred (d-O-f-) 


Xred (d-O-f-) 


Xred (d-O.f.) 


01 20053A 


5728 


20920 


621 


3-22 


20-200 


0.63 (62) 


0.59 (61) 




02 20053A 


22400 


18600 


521 


3-22 


20-200 


1.46 (62) 


0.80 (60) 


0.85 (61) 


03 20053A 


9163 


19760 


651 


3-22 


20-200 


1.07 (62) 


0.62 (61) 


04 20053A 


11040 


20550 


736 


3-22 


20-200 


2.53 (62) 


0.55 (60) 


0.55 (61) 


05 20053A 


8576 


21060 


848 


3-22 


20-200 


2.84 (62) 


1.14 (60) 


1.14 (61) 


06 20053A 


13472 


21740 


1015 


3-22 


20-200 


5.62 (62) 


0.64 (60) 


0.63 (61) 


07 20053B 


5698 


27960 


1675 


3-22 


24-200 


2.12 (59) 


1.21 (58) 


1.27 (57) 


08 20053B 


9680 


25820 


1054 


3-22 


20-200 


1.89 (62) 


0.99 (60) 


1.04 (60) 


09 20053B 


11472 


21870 


642 


3-22 


20-200 


1.03 (62) 


0.93 (60) 




10 20053B 


13920 


19320 


512 


3-22 


20-200 


1.35 (62) 


0.54 (60) 




11 20053B 


18720 


19930 


613 


3-22 


20-200 


1.65 (62) 


0.68 (60) 


0.55 (61) 


12 30036 


4208 


28590 


1840 


3-22 


24-200 


2.39 (59) 


1.41 (58) 


1.40 (57) 


13 30036 


4336 


26930 


1271 


3-22 


22-200 


1.91 (61) 


1.32 (60) 


1.42 (60) 


14 30036 


7072 


23480 


799 


3-22 


20-200 


0.85 (62) 


0.86 (61) 


15 30036 


4736 


19130 


492 


3-22 


20-200 


0.60 (62) 


0.59 (61) 




16 30036 


10320 


20000 


633 


3-22 


20-200 


2.36 (62) 


0.65 (60) 


0.67 (61) 


17 30035A 


6880 


23680 


84 


3-22 


24-200 


1.11 (59) 


1.07 (58) 




18 30035A 


6608 


19260 


43 


3-22 


24-200 


0.69 (59) 


0.68 (58) 




19 30035A 


8704 


17380 


39 


3-22 


24-200 


1.17 (59) 


0.87 (58) 




20 30035A 


15456 


18020 


50 


3-22 


24-200 


2.62 (59) 


0.71 (57) 


0.72 (58) 


21 30035A 


10064 


20040 


92 


3-22 


24-200 


2.84 (59) 


0.50 (57) 


0.53 (58) 


22 30035B 


3344 


24310 


100 


3-22 


24-200 


1.03 (59) 


1.04 (58) 




23 30035B 


7344 


16770 


39 


3-22 


22-200 


1.06 (61) 


1.08 (60) 




24 30035B 


4976 


16520 


43 


3-22 


22-200 


1.12 (61) 


0.99 (60) 




25 30035B 


16144 


18900 


76 


3-22 


24-200 


2.56 (59) 


0.49 (57) 


0.44 (58) 


26 40020 


3401 


23380 


708 


3-22 


24-200 


1.86 (59) 


0.83 (57) 


0.86 (58) 


27 40020 


3642 


20360 


544 


3-22 


20-200 


1.26 (62) 


0.81 (60) 


28 40020 


18493 


15810 


300 


3-22 


20-200 


1.14 (62) 


0.69 (60) 




29 40020 


11967 


14920 


286 


3-22 


20-200 


1.29 (62) 


1.18 (61) 




30 40020 


5735 


16100 


383 


3-22 


20-200 


0.85 (62) 


0.67 (61) 




31 40020 


7537 


17100 


470 


3-22 


20-200 


3.59 (62) 


0.70 (60) 


0.66 (61) 


32 70014 


822 


24440 


948 


3-22 


22-200 


1.16 (58) 


0.91 (56) 




33 70014 


5668 


18740 


368 


3-22 


20-200 


1.24 (58) 


0.72 (56) 




34 70015 


976 


15250 


280 


3-22 


20-200 


0.63 (58) 


0.63 (57) 




35 70015 


1248 


16300 


363 


3-22 


20-200 


0.75 (58) 


0.69 (57) 




36 70014 


3804 


17570 


447 


3-22 


22-200 


1.42 (58) 


0.73 (56) 


0.72 (57) 


37 80021 


5344 


13030 


12.8 


3-22 


28-200 


0.73 (53) 


0.59 (52) 




38 80021 


6752 


9900 


2.0 


3-22 


28-200 


0.60 (53) 


0.46 (52) 




39 80021 


17312 


8828 


7.6 


3-22 


28-200 


0.54 (53) 


0.51 (52) 




40 80021 


14624 


9071 


27.0 


3-22 


28-200 


0.64 (53) 


0.57 (52) 




41 80021 


29504 


8990 


11.7 


3-22 


28-200 


1.27 (53) 


0.78 (52) 




42 80021 


30384 


10320 


6.6 


3-22 


28-200 


1.37 (53) 


0.67 (53) 


0.70 (52) 


43 80021 


18656 


10730 


11.0 


3-22 


28-200 


1.57 (53) 


0.74 (53) 


0.69 (52) 



Note. — Wc report in this tabic, from the left to the right column, the chosen identifying number of the CDs selected spectra (as shown 
in Figure together with the associated observational dataset, its exposure times, the relative background subtracted countrates for the PCA 
and HEXTE, the energy band considered in the fits and the reduced values obtained from the fits, adopting the DBBTT, the DBBTTPEG and the 
DBBEQPAIR model; the degrees of freedom (dof) for each fit are shown in parenthesis. 



getting worse and worse from the bottom of the FB 
to the upper NB. From the residuals with respect to 
this model, it is clearly seen an excess of hard emission 
at energies above 30 keV. We report in Figure [9] four 
representative deconvolved spectra for the case in which 
the DBBTT model does not give an adequate description 
to the broad-band spectral behavior, choosing three 
spectra from the HB/NB zone and one for the topFB 
zone: spectrum 06 (at the top of the left track of dataset 
20035 A), spectrum 16 (at the top of the left track of 
dataset 30036), spectrum 31 (at the top of the left track 
of dataset 40020) and spectrum 07 (at the top of the 
right track of dataset 20035B). 

In the following we will group the spectra according to 
the relative position in their CDs, in: topFB spectra 
(spectra 7, 8, 12, 13, 17, 22, 26, 32, 37), FB spectra (1, 
9, 14, 18, 27, 28, 33, 38), NB spectra (2, 3, 10, 15, 19, 20, 
23, 24, 29, 30, 34, 35, 39, 40) and HB spectra (4, 5, 6, 



11, 16, 21, 25, 31, 36, 41, 42, 43). We caution the reader 
that this distinction is based only on the results of our 
spectral analysis for spectra taken in particular zones 
of the Z-track, without any reference to the temporal 
behavior of the source. Although it is clear, looking just 
at the lightcurves, to distinguish within one dataset the 
FB from the NB, the distinction between the NB and 
the HB is rather subtle, as the two track mostly overlap. 
For this reason a HB spectrum, for this analysis, 
is to be considered only as a simple label for a spec- 
trum extracted from the top part of left track of each CD. 

6. THE HOT COMPTONIZATION MODEL 

6.1. Energetics 

In the following we will discuss the fit results based on 
the adopted model DBBTTPEG. Although in some spectra 
this model over-parameterizes the fit, the spectral behav- 
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ior of the source and our conclusions would not change, if 
we used the best-fitting values of the DBBTT modelization 
(we show in Table [3] and d] the best-fit values of the spec- 
tral parameter, together with associated errors quoted 
at 90% confidence level). In some cases we had to freeze 
some parameters in order to make the fit physically con- 
sistent and avoid solutions without physical meaning. In 
spectra 7, 8, 12, 13, 17, 37, 38 the electron temperature 
had to be fixed to 10 keV (a value consistent with what 
we found in other CDs for the same CD zone; if free the 
fit would give a best-fit with kTg > kTe which is clearly 
wrong), the photon index of the PEGPWRLW component 
in spectra where this component is not significantly de- 
tected, or when its value is unconstrained by the fit, is 
frozen to a reference value of 2; the DISKBB parameters 
of spectra 42 and 43 had to be fixed to the best- values as 
obtained by the best-fit model DBBTT (if free, this com- 
ponent vanishes, while the COMPTT and PEGPWRLW have 
comparable flux values). 

In this spectral decomposition we calculated the total 
luminosity as follows: 

= A^D'^{Ftt + FpEG + a^NoBTU) (1) 

where Ftt is the COMPTT flux and Fpeg is the 
PEGPWRLW flux extrapolated in the 0.1-200 keV energy 
range; g — 1/(/^k^), where / is t he spectral hardening 
factor (jShimura fc TakaharalllQQSh . is assumed equal to 
1.7, and k ~ 0.4 corrects for the fact that the radius at 
which is reached the maxim um temperature is grea ter 
than the inner disk distance (jVierdayanti et ani2006f l. a 
is a constant dependent on the distance of the source, 
which we assumed to be 2.8 kpc, and the inclination of 
the disk, which we assumed 45°. In Figure[l]we show the 
contribution of the disk and of the Comptonized compo- 
nent to the total flux. The luminosities on both axis are 
in Eddington units for a 1.4 Mq. 

The source span a luminosity range between 0.9 IjEdd 
and 2.4 IjEdd, thus being in a super-Eddington state al- 
most independently of the position of the source on the 
CD. We note that the only way to drastically reduce this 
range is to suppose a much massive compact object (the 
possibility of a massive NS h osted in the Scorpius X- 
1 system is also discussed in IStella fc Vietril ()1999[ ) in 
the framework of the relativistic precession model used 
to explain the frequencies of quasi-periodic oscillations 
in NSs and BH systems) given that both the distance 
and the inclination angle are known with good accuracy 
(10% and 13% relative error, respectively). 
The disk flux and the hard flux correlate as expected 
with the total flux, although the disk emission presents 
a more pronounced linear relation, while the hard flux 
seems to saturate as the total flux reaches the highest 
values. If we assume a linear dependency of the two con- 
tributions versus the total emission we derive a constant 
of proportionality equals to 0.56 for the disk flux and 0.44 
for the Comptonized flux. When we look at the ratio of 
the two contributions along the Z-track (see Figure [2l 
a 5% error on this ratio is assumed, after having tested 
the corresponding hard and soft flux values for different 
parameter values of the spectral components inside the 
uncertainties errorbars as derived from the fits), there is 
a general trend to have a higher hard/soft ratio as the 
source moves from higher to lower accretion states; this 
is in agreement with our expectation that, at lower ac- 



cretion rate the disk could be truncated and the fraction 
of the power dissipated in the boundary layer should cor- 
respondently increase. Spectra 7, 12, 13 and 32 do not 
follow this trend; these spectra belong to the topFB spec- 
tra, and in these cases, we are possibly, underestimating 
the disk contribution, or overestimating the hard fiux. 

6.2. Disk emission and thermal Comptonization 

We present in Figure[3]the dependency of the inner disk 
radius and of the inner disk temperature, plotted against 
the hard/soft luminosity ratio. It is shown, both in Ta- 
ble [3l [4] and in Figure [H the apparent radius, derived 
without considering important correction factors, due to 
the color/effective temperature ratio, relativistic effects, 
and the non-zero torque boundary conditions; moreover, 
it should be noted that these corrective factors can also 
vary, depending on the accretion state. The sum of all 
these effects could be a drastical rescale of the measured 
inner radii up to al most an order of magnitude higher 
(jMerloni et al.l [20001 ). but the general trend that corre- 
lates this parameter with the accretion state of the source 
should be qualitatively preserved. 

The apparent inner disk radius, for spectra in the topFB 
has the lowest values, < 10 km, while the average value 
for NB spectra and HB spectra is considerably higher 
(average value: ~ 18 km for the HB and NB, ^ 13.5 
km for the FB and ~ 10.2 km for the topFB). The in- 
ner disk temperature, on the other hand, presents the 
highest values, as one could expect, on the FB, while the 
temperature diminishes with the accretion rate. A small 
subset of FB spectra falls out of the trend, but as pointed 
out earlier, in this case, the exact determination of the 
relative contributions of the two components could suf- 
fer of systematic uncertainties linked to the particularly 
luminous state of the source. On the FB the disk tem- 
perature has an average value of 2.12 keV, while on the 
lower part of the FB sinks to 1.73 keV. There is not any 
appreciable difference on temperature for spectra taken 
on the NB and on the HB with an average value of 1.43 
keV. 

For the hard, thermal Comptonized, component the 
accurate determination of all the spectral parameters 
was not always possible. The soft-seed photon tempera- 
ture presents, for each spectra analyzed, a substantially 
higher value with respect to the disk temperature at the 
inner radius. We hence propose that the soft photons are 
originated from the boundary layer/NS surface, partially 
thermalized with the softer photons of inner part of the 
accretion disk around the NS. Moreover, calculating the 
radius of the soft-seed photons emitting region from: 



R^ = 3xl0UJ-^/ikTof, (2) 

where d is the distance to the source in kpc, fboi is the 
COMPTT bolometric flux and y (see Eq. [3]), the Compton 
parameter, we derived for our sample values in the 3-6 
km range, that clearly indicate a rather small emitting 
region, thus supporting our identification. 

The CD correlated changing of the thermal temper- 
atures is compatible with a scenario in which the posi- 
tion of the source on the CD is determined by the in- 
stantaneous M, and higher accretion rates correspond to 
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Fig. 2. — Left panel: soft (disk) luminosity for the extracted spectra vs. total (0.1-200 keV) luminosity. Right panel: hard 
(Comptonized) luminosity for the extracted spectra vs. total (0.1-200 keV) luminosity. All luminosities are in Eddington units 
for a 1.4 Mq NS. 



12 13 1i 



Fig. 3. — Ratio of the hard/soft luminosities for all the CD extracted spectra. In the x-axis the number refers to the spectrum 
as indicted in Table (2] We assumed a 5% error for ratio of the two fluxes. 





inner disk radius vs. the hard/soft ratio. Right panel: temperature at the inner disk radius vs. the 



Fig. 4. — Left panel: 
hard/soft ratio. 

higher seed-photons temperatures, resulting in a hotter 
radiation field. The soft seed-photon temperature has a 
2.2 keV, 2.4 keV, 2.9 keV average values for spectra on 
the HB/NB, FB and topFB respectfully. When plotted 
against each other, the disk temperature and the soft 
seed-photon temperature (see Figure [SJ left panel) follow 
each other quite closely, being these two parameters the 
main driving physical quantities related to the accretion 
state of the source. The correlation between the CD re- 



solved spectra and the two temperatures shows, although 
in a qualitative representation, that the link between the 
spectral evolution of the source and its accretion state is 
well motivated. 

The other two spectral parameters, that describe the 
high-energy curvature of the Comptonized component, 
the electron temperature of the Comptonizing plasma 
(kTe) and the optical depth of the cloud r, are not al- 
ways well constrained by the fit. As it is evident from 
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Hard luminosity (Edd. units) 



Fig. 5. — Left panel: soft seed-photon temperature kTo vs. the inner disk temperature kToB- Right panel: soft seed-photon 
temperature kTo vs. hard (COMPTT+PEGPWRLW) luminosity. 




Fig. 6. — Left panel: optical depth vs. electron temperature (all spectra). Right panel: optical depth vs. electron temperature 
(subset of spectra on the HB/NB). Errorbars show 1.0 a errors. 



Figure [6l for spectra at high accretion rate we obtained 
high values for the electron temperature (kTg > 10 keV) 
and low values for the optical depth (r < 1). As the 
sources resides on zones of lower M , the Comptonizing 
cloud substantially thickens, while the electron temper- 
ature correspondently diminishes. Plotting the subset of 
the spectra in the NB and HB helps to better visual- 
ize this tendency (Figure [6l right panel). It is to be 
noted that the determination of these two parameters is 
partially related with the hard power-law component; al- 
though there are spectra where the PEGPWRLW component 
is not strictly necessary for a good fit, its introduction 
slightly shifts the the thermal Comptonized component 
to lower energies, thus giving a noticeable reduction in 
the uncertainties of the Compton curvature. 
In all the spectra it is required a Gaussian line at ener- 
gies ^ 6.4 keV. It mostly appear rather broad, with a line 
width that often reaches our constraint of 0.8 keV. The 
equivalent widths associated with the line are in 80-200 
eV range, and are to be taken as rather indicative given 
the low energy resolution of the PCA and the presence 
of systematic residuals in this energy range. 

6.3. The hard tail behavior 

The presence of a hard X-ray excess in the extracted 
spectra of Scorpius X-1 is mostly evident in 18 spectra 
from a total of 43 CD selected spectra. There is little 
evidence of the presence of hard X-ray excess for spec- 
tra lying in the bottom part of the FB and of the NB, 
namely for spectra taken near the apex that connects the 
left and right track that compose the V pattern. 
In all the examined CD patterns we always detect a 
power-law high-energy excess, dominant above ~ 40 keV, 
as the source resides at the top of the left track, i.e. a 
portion of the diagram that we tentatively identified with 



the HB, and which should correspond to the lowest mass 
accretion rate. This group of spectra are characterized 
by the lowest values of inner disk temperature and soft 
seed-photon temperature; the derived inner disk radii 
correspond to a disk truncated at about 10-15 Rg, while 
the Comptonizing optical cloud is substantially thicker 
with respect to any other spectrum in other zones of 
the CD; the fraction of the total power dissipated in the 
Comptonizing corona has the highest values, with a ratio 
hard/soft luminosity well above unity. 

Apart from this consistent group of HB/NB spectra, 
we also found a minor group of spectra (namely 5 spec- 
tra: 7, 8, 12, 13 and 26), for which the fit required a third 
high-energy component. These spectra belong to 3 CDs 
(2GG35B, 30036 and 40020) and are spectra taken at the 
top of the FB. Although these spectra are located in the 
same CD region, not all the CD that display a topFB 
clearly show this excess. This could be due to a statisti- 
cal reason, as in all the other CDs the HEXTE countrate 
(see Table [2]) drops of about one order of magnitude 
with respect to these 3 CDs; because the countrate is 
essentially concentrated in the 20-35 keV energy band, 
that is the range that constrains the curvature of the 
thermal Comptonized component, the energy channels 
in this energy band with a high statistic will drive the 
fit, and any small excess above the exponential tail of 
the Comptonized component will statistically be more 
accentuated. However, we cannot exclude that an oppo- 
site reasoning is true, i.e. that we are underestimating 
instrumental effects that appear only for higher count- 
rates as dead-time corrections or poor background esti- 
mate, so that these particular hard tails are artifacts of 
an incorrect modelization of the thermal Comptonization 
curvature. As we have no apparent way to better cali- 
brate our spectra, we discuss also the appearance of this 
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Fig. 7. — Left panel: hard tail luminosity (in the 20-200 keV energy range) vs. the hard/soft ratio. Right panel: the energy 
index of the PEGPWRLW component vs. its 20-200 keV luminosity (the subset of spectra for which the photon index was not 
costrained by the fit is excluded from this figure). 

component for this group of spectra. Future observations 
with other high-energy satellites, such as INTEGRAL or 
SUZAKU, will test if this particular topFB state of the 
source is also accompanied by a high energy excess or it 
is an HEXTE faked detection. 

The hard excess has no evident relation with the total lu- 
minosity of the source, but this could be expected since 
the different zones of the CD track are not related with 
the total X-ray luminosity. Using the pegged power-law, 
we directly derived from the normalization value of this 
component, the flux in the 20.0-200.0 keV energy band. 
We show in Figure [51 left panel, the hard tails lumi- 
nosity in the 20-200 keV range vs. the hard/soft ratio. 
From the plot it can be clearly seen that the hardening 
of the spectrum is reflected in the increasing luminosity 
of the hard-tail component; the CD resolved spectra are, 
correspondently, disposed on this plot: the FB spectra 
reside on the bottom left of the diagram with hard fluxes 
< 2 X 10~^° erg s~^, NB spectra are harder and half of 
them show significant fluxes 4 x lO""'^" erg s~^), while 
the HB spectra, in the right part of the plot, show the 
strongest hard X-ray emission (in the range 6-12 xlO"^*^ 
erg s~^). 

As pointed out in the previous section, we found for the 
topFB sample, spectra with no detectable hard X-ray 
emission, that occupy the left bottom part of this figure 
and that smoothly join the correlation between hard/soft 
ratio and hard flux and spectra that do not follow this 
trend and for which there is both in the CDMPTT and 
in the PEGPWRLW component a hard flux higher than we 
expected. From the figure is also evident that one of the 
most luminous hard tail is found in a particular spec- 
trum on the topFB (namely spectrum 26 of CD 40020); 
the luminosity of this component is related to the to- 
tal broad-band high luminous state of the source, that 
reaches in this case the highest value of total luminosity 
(2.6 hEdd) observed in our sample. 

On the other hand, the values of the photon index present 
a well defined bimodal distribution according to the po- 
sition of the source on the CD. We find quite flat power- 
laws, with index values less that unity for all the FB and 
topFB spectra, while for the other group of spectra, val- 
ues generally range between 1.5 and 2 (Figure [SI right 
panel). 



6.4. The hybrid Comptonization model 

Models of hybrid Comptonization have so far been 
mostly adopted to explain state transitio ns both in black- 
hole candidate systems (as in Cyg X-1. iGierlinski et al.l 



[Tool GRS 1915-H05, IZdzia.rski et al.l [2OOII: or GX 

339-)-4, Wardzinski et al. '200^ and in NS system (GX 
17+2, .FarineUi et al..,2005,) . 

We tried to model the spectra for which we detected a 
hard tail, making the hypothesis that both the hard tail 
and the thermal Comptonized component were related 
to each other. In this way, all the broad-band evolution 
of the source is covered in a self-consistent way by only 
two spectral components: the thermal soft disk emission 
and the hard hybrid Comptonized component. In table 
[5] we present the results of our fits for the group of 18 
spectra, distinguishing between the HB/NB group (top 
part of the table), for which we used the DISKBB compo- 
nent as in the case of DBBTTPEG modelization, and the 
topFB spectra, for which we used the DISKPN compo- 
nent, under the assumption that for this state the disk 
is not truncated, and reaches the surface of the NS (we 
assumed that this happens at a distance of 6 Rg). 

As in the case of the previously adopted modelization, 
the thermal disk temperature and the soft-seed photon 
temperature are the main driving physical parameters 
that determine the changing of the spectral state of the 
source. We show in Figure [5|(left panel), the correlation 
between the two temperatures: the topFB spectra 
occupy the top right part of the plot, with an average 
kTBB value of 2.73 keV and inner disk temperature 
kToB of 1.97 keV, while the set of the HB-NB spectra 
are disposed along a linear trend in the 1.2-1.8 keV 
energy range for the disk temperature and 1.9-2.5 keV 
range for the soft seed-photon temperature. The plot 
closely follows the plot of Figure [51 right panel. 
The hard compactness 4 results for all the examined 
spectra only a small fraction of the soft compactness 
Is, with values of hard/soft ratio Ih/h below 0.2. The 
HB/NB spectra span almost an order of magnitude in 
the Ih range, while the topFB spectra are all grouped in 
a more narrow range of values, around 4 — 0.6-0.8 (see 
Figure [51 right panel). 

The two groups of spectra significantly differ both in 
the T and Inth/lh values (it is to be noted, however, 
that the optical depth reported from the fits has not 
the same physical meaning of the classical Thomson 
optical depth of the Comptonizing cloud, as it more 
properly refers to the optical depth associated with the 
background photon radiation field): HB spectra present 
generally higher values of optical depth (with an average 
value of Tp ~ 2.2) and a hybrid electron distribution 
(with 0.3 < Inth/lh < 0.8, and an average value of 0.55), 
while the topFB spectra have considerably lower values 
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Fig. 8. — Left panel: soft seed-photon temperature kTss vs. inner disk temperature kTmsK (which corresponds to the 
kTi5_B value of the DISKBB component for the HB/NB spectra and to the kT,nax value of the DISKPN component for the topFB 
spectra) for spectra fitted using the DBBEQPAIR modelization. Right panel: the soft seed-photon temperature is plotted against 
the hard compactness parameter. 



oi Tp < 1 and essentially non-thermal spectra, with 
Inth — 1- The slope of the high-energy power-law of the 
hybrid electron distribution, Tinj is consistent with a 
value of 2, for the HB-NB spectra, while for the topFB 
spectra its value is considerably lower (between and 1). 



7. CONCLUSIONS 

We outline our main conclusions from an extensive 
analysis of Scorpius X-1 RXTE observations in the 3- 
200 keV energy range: we observe a spectral evolution 
of the source that is clearly dependent on the position 
of the source along its CD track; the spectral decompo- 
sition consists of four different spectral components that 
we interpret, from the softest to the highest X-rays re- 
spectively, as follows: a soft thermal component from an 
optically thick accretion disk, a thermal Comptonization 
from the boundary layer, a reflection component which 
mainly manifests in a broad Gaussian line at 6.4 keV and 
a variable hard excess above 30 keV, mostly present at 
low M. 

We have shown that the CD correlated variations of the 
two thermal temperatures (disk temperature and soft 
seed-photon temperature), the trend in the hard/soft ra- 
tio, the Compton thickening of the Comptonizing cloud 
and the apparency of significant hard X-ray emission on 
the HB/NB are the main spectral characteristics of the 
source. 

We paid particular attention to the modelization of the 
hardest X-ray component, which is of primary impor- 
tance for the understanding of all the continuum emis- 
sion, as it has a major impact on the determination of all 
the other spectral parameters. To fit this hard excess we 
modelled our spectra using a phenomenological compo- 
nent, i.e. a power-law with a low energy exponential cut- 
off at the seed-photon temperature and a self-consistent 
physical model, a hybrid thermal/non-thermal Comp- 
tonization code. In both cases we found an adequate 
description of the spectra for each source state, and the 
two modelizations are statistically equivalent. 
The use of simple power-laws to fit th e hard tails in the Z- 
class sources, but also in some atolls (jTarana et al.l[2007l : 
iFiocchi et al.l 12006). has been largely used in the past, 
given the lack of broad-band coverage up to the MeV en- 
ergies, and a poor understanding of the physical origins 
of this component. C ontrary to what previously reported 
p'Amico et al.ll200lD . the presence of the hard tail is re- 
lated to the broad-band spectral evolution of the source 



as inferred from the source position on the Z track of 
its CD: the presence and the values of the photon in- 
dexes in Scorpius X-1 are consistent with the results 
from fits to spect ra on the same zone of the CD for the 
other Z sources (jAsai et all Il994t iDi Salvo et "all l2000l . 
[200ll[200il : at the same time, the flux contribution to 
the total energy output of the sources is of the same or- 
der and is anticorrelated to the inferred mass accretion 
rate. This is clearly shown in the case of the CD 20035A, 
where the hard tail flux monotonically increases as the 
source moves from the bottom to the top of its NB/HB. 
The only difference, between Scorpius X-1 and the other 
similar sources, is constituted by rather flat hard tails 
that we found when the source is at the top of the FB 
and which has never been observed in other NSs sources; 
however, while for every CD that we obtained, we sys- 
tematically found a hard tail as the source was at the top 
of left V track, the detection on the FB is presumedly 
dependent on the HEXTE statistics. Past surveys with 
the use of scintillation counters on-board on balloons re- 
ported the unusual flattening of the spectrum of Scorpius 
X-1 above 40 k cV ( Agrawal ct al. 1971), as well as more 
recent surveys (|Manchandall200Q) : it would be interest- 
ing to definitively assess the existence of such compo- 
nent; from our analysis, we can only conjecture that this 
component, could be a signature, on the FB, of the tres- 
passing of the Eddington limit on the interface between 
the accretion disk and the board of the boundary layer, 
followed by a violent expulsion of part of the accreting 
matter. 

Hard tails on the HB could have a thermal origin, as 
in the case of the thermal Comptonized component that 
dominates the spectrum at lower energies (less than 20 
keV); in this case the power-law, as long as the optical 
depth is low and the electron temperature is not too high, 
is a good approximation of the Comptonized spectrum, 
whose y parameter 



y - 



X (t + t^) 



(3) 



is related to the photon index of the power-law by the 
following relation: 



a 




(4) 



In the case of Scorpius X-1 spectra on the HB/NB, we 
derived photon indexes in the 1.8-2.2 range, that imply a 
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Fig. 9. — Deconvolved spectra, in the E*f(E) representation, and contributions per channel to the x*^ for four spectra of our sample; first 
column: spectrum 06, second column: spectrum 16, third column: spectrum 31, fourth column: spectrumOT. Deconvolved spectra (upper panels) 
and contributions to the (lower panels) from the best-fit using the DBBTT, DBBTTPEG and DBBEQPAIR models for the first, second and third row 
respectively. In the spectral decomposition the green line component shows the disk emission, the blue line shows the Comptonized emission, red 
line the Gaussian line, the thin blue line the PEGPWRLW component. 
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y Compton parameter in the 0.6-1.1 range; 'at — 1 thick 
plasma would then require a ~ 50 keV thermal plasma, 
while a T = 2, on the contrary, a ~ 15 keV electron tem- 
perature. We did not find any evidence of high-energy 
cut-ofi^, obtaining only lower limits in the 50-100 kcV 
energy range, and this constraints the optical depth to 
be less than one while the electron temperature would 
be greater than 50 keV. The formation of hot zones, or 
blobs, of very hot plasma, possibly powered by episodes 
of magnetic reconnect ion above the d isk, constitutes a 
plausible scenario ((Haardt et al.l 1199^1 ): this model has 
been also prop osed in the case of t he hard/low states of 
BH candidates (|Malzac et al.l[200^ , which share with the 
bright Z-sources systems, from a spectral point of view, 
similar values of the hard photon indexes and the low ac- 
cretion state, while they naturally differ from the strong 
boundary layer emission present in the NSs sources and 
absent in the BH case. 

A thermal origin, on the other hand, for the topFB hard 
tails seems to be excluded, as it would require unrealistic 
large values of the electron temperature. 
Our analysis has shown that this scenario is statistically 
equivalent to a hybrid Comptonization model, where the 
non-thermal fraction of the power injected in the electron 
heating counts up to half of the total injected power. The 
non-thermal fraction is significantly highe r with respect 
to the value found bv lFarinelli et al.l (j2005l ) in the case of 
GX 17-1-2. The non-thermal fraction for the topFB spec- 
tra must be distinguished from the non-thermal spectra 
of the BH sources, as in the case of Scorpius X-1, the 
hard luminosity is a small fraction of the soft luminosity, 
that dominates the 10-30 keV energy spectrum. Most 
of the soft photons do not Compton interact with the 
electron cloud, because of the low optical depth values. 
Another possible physical mechanism to explain the hard 
X-ray variability is the bulk motion Comptonization (e.g. 
I^saltis 2001), present in systems with a high velocity ra- 
dial accretion flows. This radiation mechanism, used to 
explain the hard power-law component in BH systems 
spectra in their soft/high states, is, however, not able to 
account for the the absence of any cut-off up to energies 
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of 0.5 MeV and for the hard photon ind exes observed 
(r < 2.0) (|Niedzwiecki fc Zdziarskll2006f ) . In the case 
of bright NS systems, the high radiation pressure of the 
inner boundary layer constitutes a strong barrier for any 
incoming convergent bulk motion flow. On the other 
side, because jets have been observed in the radio from 
Scorpius X-1, bulk motion Comptonization inside the jet 
can be an important production mechanisms also of hard 
X-ray radiation. Theoretical spectra and energetic con- 
tribution according to the relevant physical parameters 
involved in the case of a strict coupling with inflow (i.e. 
accretion to the compact object through the formation of 
an accretion disk) and outflow (i.e. jets), are a promising 
way to cover in a self-consistent way all the phases of ac- 
cretion. Details about the geometry of the scattering me- 
dia, amount of reflection on the cold disk and dependence 
from the accretion state of the source are yet to be fully 
explored. Moreover, in this case synchrotron emission 
of soft photons by a beamed population of relativistic 
electrons becomes a competitive source of photons with 
respect to the thermal soft photon emission of the bound- 
ary layer and accretion disk. An attempt to explicitly 
compute a jet spectrum, fr om radio to ha r d X-r ays, has 
been recently proposed bv iMarkoff et al.l (|2005f) . where 
the jet base subsumes the role of the static Comptonizing 
corona; spectral fits in the case of BH systems (namely 
Cyg X-1 and GX 339-4) in hard states are consistent with 
this scenario, but BH soft states and NS systems spectra 
need yet to be tested in order to understand the limits 
of validity of the jet model. 

Questions to be further addressed in future observations 
are: the exact shape and contribution of the soft compo- 
nent below the 3 keV range, the investigation with good 
spectral resolution in the 6-10 keV range of the reflection 
features, the extension up to 0.5 MeV of the spectral cov- 
erage in order to constrain the physical parameters char- 
acterizing the non-thermal electron distribution of the 
Comptonizing plasma. The latter point is well within 
the capabilities of the recently launched SUZAKU satel- 
lite, so that future observations in this direction can be 
a stringent test to our conclusions. 
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-'-"•''-1.0 

i6.5±?;^ 


6.40+"- 

6.40+O-is 
6.44i;!;l? 

6.49^0.17 
6 52+01^ 


0+o.s 

0+0.4 

r, + 0.3 
'J. ■'-(). 3 

o.sli;:^ 

0.57l°J 


0.08llil« 
0.09in^ 

o.i4i!i:;;i 

0.20lO:»f, 

g 17+0. 08 
'-'•-'-'-0.06 


2.92i;!:iii 
2.801;!:- 
2. .34+11"/.; 

2 l2+"'o-^ 
9 09+0.05 
^•''^-O.OS 


10 

10 

q+2 

•'•'^-0.4 

4 50+2 


0.021;;:- 

o.oil'^-i* 

5I3 

6.311;; 

3.89lg;?5 


"••^-0.3 

5+0-5 
'-'-3-0.5 

'J-^-0.2 

-■-••^-0.5 

1 4+0.5 


-30+180 

500+"0 
^'^'^-120 

26oi?,;; 

390l|'(; 
390+?° 




10.8 
7.76 

7.53 
8.3 
11.4 


12 
13 
14 
15 
16 


9 fi7+0.U9 
-0.09 
9 nc+O.OS 

^■UD_Q gg 

-i_n Hzi 

l.86l°;»^ 

J--48-o.Ofi 
1 42+0.03 
-^•^^-0.06 


»-'-0.7 
Q 0+0-5 

1 n f^+o-b 

15 7+"^''^ 
-'-■J- ' -0.8 

ir.olii 


6.40+0.18 

6.40+0.21 
-i_n 1 ^ 

fi 4S+0.15 


0+0.5 
0+0.5 
0+0.5 

o.slS:^ 

0.65l°;« 


09+" "'' 
n 09+0.04 

U.U9_Q Q5 

nil +o.ue> 
'J.ll-0.04 

n 20+0 08 

'J-^'J-0.07 

o.2ol°:°? 


9 QO+U.U2 
^•='•-'-0.07 
2 85+0 03 

-i_n n"^ 
9 c;fi+o.oo 

^••J'J-0.06 

9 90+0.05 

^■^•-'-0.04 

9 01 +0.03 

^•■^^-0.03 


10 
10 

-1-1 

9+f 

4.9+1^ 

t; 4+0.7 
^^•^-0.5 


0.01+0.16 

0.01+0." 

-LA 

0.61^3 

3.211.4 
9 7+1.5 
^•'-0.1 


7+0.6 
^•'-0.6 
0+0.3 

2 

1.201°;*' 


560ll«° 

48011?° 
-LI nn 
^'-'-50 
1001?°° 

6001^° 




12.1 
10.1 
7.59 
7.15 
8.54 


17 
18 
19 
20 
21 


9 -, c:+0.U3 

"^-0.06 
1 cq+0.05 
-•-•"J •-'-0.07 
1 Q4+0.13 
■'■•'^^-0.12 
, 90+0.22 
^•^•^-0.12 


ii.7±i:l! 

15.511- 

i8.6l^:» 

2411 

soil 


6.40+0.26 

f, 40+0. 16 
O.4O-0.0.8 
fi 51+0.15 

fi 42+0.25 
"•^^-0.02 
f. 49+0.25 
"•^^-0.02 


0+0.5 
0+0.3 

5+0.3 

0-8_o.4 
0-8_o.4 


o.i8l{l:?« 

24+0." 
'J-^^-0.08 

26+0-" 
^•^'^-0.08 

41+0-" 

"•^^-0.19 

^2+0-l'5 
"•^^-0.29 


9 09+0. 05 
^•°^-0.07 
2 40+0." 

^.■''^-0.03 
9 -1-1+0.16 
^.l-'--0.17 
1 QP;+0.37 


10 

3.6+? 

4 3+^ 
7+0.7 
•^.'-0.3 
9+0.5 
•^•^-0.2 


0.01+0.19 

5+^ 
^-3 

4.8+1-* 
6.2+1-'' 
9.9+1-^ 


2 
2 
2 

-, 7+0.5 
1- '-0.4 
9 9+0.4 


0+100 

190111° 
390111° 
7201^^0 

1540111° 




11.1 
10.4 
10.5 
13.8 
21.8 



TABLE 4 

FIT RESULTS: THERMAL COMPTONIZATION MODEL 



Parameter kT^s Res Egaus dgaua t^gaua To Te r Index Fpeg Ftt 

Units keV km keV keV ph cm~^ keV keV erg cm~^ erg cm~^ 



22 2.20t°o\i lltlo 6.40+03 0+»5 0.20«-J2 2.90l:lil lltf 0.14^23^3 2 lOOt^OO 14.0 

23 1.56+_rit 17.2lli 6A3+_l:ll 0.751";^ 0.45«;}| 2.30i»:i^ 4.otJ-i ^-O^li 2 1201^ 9.8 

24 l-30tH? 24tf 6.49tE;-20 0.68t°-i2 0.30lO-}t 2.00t»i 3.7«I 6.3i°i 2 SOOl^go 13.2 
25 1.48t°-[i^ 19.4±^-^ 6.54lg'j^ O.eslgj^ O-Set"'}^ 2.42±»-°3 5.2^^-^ 3.5^11 l.elg'^ 940i;i° 16^ 

Te TT^, UEI^, ^W^, ^^^o ^^^^s ^7 ^^^^ ^ 

27 2.0210-08 11.7lE;| 6.5410^0 0.48t0-3 0.29«-}6 2.821H6 OS^lt 1.8l°? 540111° 10.7 

28 1.60io;« 14.6ti;o 6.50t»;« o.6lt°i« 0.29«;j? 2.30tH^ l^lig O.Oel"-^ 0.5lj;° 250l|° 8.2 

29 1-441^4 17.1110 6.50l°;3 0.68lE!-^2 0.26tO-«? 2.25in3 ^-^-ol ^-^-l 2 2401^20 g.S 

30 l-45iS-g^ 17.0tl-6 6.49tE!-2i 0.80_o 3 0.28tO-;i ^.SSt^'ol ^■'^-17 2.4t? ? 2 270l2;50 9.6 
31 1.46tg;;^ 18.5llo 6.521°;^^ 0.5tg;3 o.25tg;;? 2.34tg;^^ 3.71^;^ 6.0l°;g 1.5l[i:^ llSOtl^g ^ 

32 2.09t^-g| 10.7+^!; 6.40+0-29 oF^^I 0.22t^-;» 2.84+^;-^^ 10 0.1lt^-*o 0:5^^^^ 600^1™ leT 

33 l-79tHI 12.711-1 6.40+0-24 0.41°-^ 0.23^-1^ 2.5tli "^-^tH 3.91^7 0.81°-^ 300±i°° 8.7 

34 l-47t°-E;^ 18.5tl-' 6.40+0-27 0.68l°-i2 0.20tO-li 2.2t1i 3.4^^^^ 7^1 2 300ll80 g.g 

35 l-40j:0-[]? 19-3t2-? 6.40+0-24 0.62l0-f 0.22tO-J3 2.3lin6 S.ltfj 3.01^-0 2 lOOl^oo 10.3 
36 1.46t0-[i^ 16.8±1-^ 6.40+''-3 0.8_o 4 0.29^0-°' 2.381^8 4.llg-^ 4.6^^-^ 1.9±i;-^ 450t^o lOJ 

37 2loT^^| ll.ltg-^ 6.47tg-^? OeT^^I 0:28T^^| 2.83tg-;;i 10 O.Olto* 2 TooT^^ ei? 

38 l-7ltH9 I3.9l?1l 6.40+0-22 o.8_o 2 0.38t012 2.4410-" 10 0.4t0-3g 2 2501^00 7.4 

39 l-45tH6 IS-Sill 6.40+0-20 0.8_o 2 0.40l0 02 2.23t0-J0 n+f o.elj-^ 2 100^00 8.6 

40 1.42to:o6 19,6+2:0 6.4210;^^ 0.8_o:2 0.35to:io 2.22lo;io lotf O.StH 2 901^° 9.7 

41 1.5110-06 17,2+1.4 6.44+0:23 o.8_o.2 0.33tO;0^ 2.34t0:0| ^tl' 0.4lgi 2 170l|0 9.6 

42 1.47 19 6.47l|l-gf 0.8_o 4 0.28^1]-?^ 2.36lg-l!l 4.6^1-^ 3.4^3 ^ 2.5lg-^ 400l^go ii,g 
43 L36 22 6.40+0-24 0.8_o ^4 0.40^0-13 2.25t°:01 4.5t^-j 4.210;^ 2.010;^ 650^^0;, 14.6 



Note. — Parameters values and associated errors for all the CDs resolved spectra adopting the DBBTT or the DBBTTPEG model. The column FpEG gives the power-law flux in units 
of 10~^2 erg cm~2 calculated in the 20-200 keV energy range. The last column Ftt gives the flux of the thermal Comptonization component in units of 10~3 erg cm^2 g— 1 in 
0.1-200 keV energy range. See Section |4] for an explanation of the parameters involved in the flts. All the uncertainties are calculated at 90% confldenco level. 



TABLE 5 

FIT RESULTS: THE HYBRID COMPTONIZATION MODEL 



Parameter 




Rrn (K) 


LineE 




^Gau3 


ih 


kTBB 


Inth/^h 


Tp 






Units 


keV 


km 


keV 


keV 


ph cm~2 g-i 




kcV 














-'-■4^_0.04 


1(5+''' 


b.5U_o 1 


'J-"°-0.23 


U.i»_o.04 


n fi>^+o.U3 


^•^'^-0.08 


00+9 

3S_g 


0-8i_o.l3 


9 


o.oy 


04 


1 O1+0.05 
^"^^ — 0.04 


19+f 


6 45+°-?,^ 


8 „ = 

"-"—0.5 


26+0-0^ 


98+0-00 


9 91 +0.06 
^-^^—0.06 


62+4 


1 7+0-7 

' -0.7 


2 


11.52 


05 


1.30+°°? 


19+9 


6.50+°-?° 


0.8 n ^ 
— 0.4 


0.22+°-f 


1.01+0-" 

^ — 0.17 


2.24+°-°5 
'^•■^^— 0.07 


63t| 


1.8+0-3 

^ — 0.7 


2 


12.33 


06 


1.44+°-°4 

— v. 12 


16+1° 

— D 


6.50+°-?° 

— U.l 


0.46+°-°? 

— U.ll 


0.12+0-04 

— U.U4 


0.88+0-14 

— 0.05 


2.46+°-°3 

— O.Oo 


73+4 


0.64+0-^„ 

— 0.09 


2 


12.46 


11 


1.48+°-°5 

^—0.14 


18+i^ 


640+0.16 


0.6+°-2 


0.25+0-" 

^ ^ — 0.07 


1.12+0-}° 

— 0.10 


2.17+°-°^ 
' —0.14 


53+11 


9 ai+0.14 


2 


13.25 


16 


1.16+°-°5 


— ID 


6.40+0-2° 


0.8 n 
" ^ — 0.3 


0.38+0-03 

— 0.07 


2.02+0-1,! 

— 0.04 


1.93+0-?f 

— 0. lo 


75+7 

— 3 


7.3+0-3 

— 0.3 


2 


11.39 


20 




23+? 


6.50+°-f 


6+°-2 
U.D_Q 2 


SI +0-2 


1.32t0.3 


2 06+°-°'' 
^•'-'0-0.07 


54+3 


4.0^0-8 


2 


14.99 


21 


1 24+°°* 




"•+"-0.06 


0-8„o.4 


4S+013 


-■--^^-0.17 


9 19+0.1 
^--'-^-0.02 


52+3 


9 + 1.6 


2 


12.31 


25 


1 24+°-°l 




fi 42+0-1** 
"•+^-0.02 


0.8_o.2 


41 +0.06 
^•+^-0.09 


-I 90+0.13 


9 oq+0.04 

^-'jy-0.03 


47i3 


9 0+0.6 


2 


11.14 


31 


1 O1+0.11 

^^•■^J^-0.02 


22+9 


f. 40+0.22 


0-8-0.2 


•^fi+0 06 


1 11+0.09 
^•^^-0.16 


9 20+0-17 
^-^'-'-0.06 


67+1 


1 Q+0.4 
1-^-0.7 


2 


16.52 


36 


1 44+0-08 




6.40+0-27 


7+°-l 
'-'■'-0.5 




0.6210-23 


9 04+0.09 


50+10 


6+0-8 
U.D_o 4 


2 


2.47 


42 


1 CI +0.06 

^•'^J^-o.os 


lll^ 


6.40+0-3 


0-8_o.5 


„+0.03 
'J-J^^-0.05 


OO+0.09 
U-J'J-0.04 


2 40+0.07 
^-^'-'-0.12 




90+0-6 
'J-^^-0.15 


2 


12.1 


43 


1 41+0.07 
-'-•+-'--0.04 


121^ 


fi 42+0-26 
"•+^-0.02 


0-8_o.3 


„ ,,c:+0.008 


64+0-11 


2 29+°-i3 
^-^^-0.06 


411* 


6+0-6 
U.D_o 4 


2 


1.4.3 


07 


1 95+o.o;i 

^■^3-0.06 


^■'J'^-0.06 


64+0.17 


0.0+0.5 


'^•"^-0.03 


0-85+0.21 


2 74+0.10 
^- '+-0.06 


100_i 


o.eot^.^^ 


Og+U-8 
"■^-0.9 


12.62 


08 




-1 17+0.05 
^•^'-0.07 


6.4+0-2 


0.0+0.5 


in+0 03 

"■J^'J-0.05 


0.89+°-°« 


9 fio+0.07 
^-"°-0.03 


100_4 


0.47+g-3 


O-SlJ.g 


8.62 


12 


1.981°;?^ 


1 nn+°-°i 


6.4+0-2 


0.0+°-4 


nq+0-03 


0.72t°iii 


9 09+0.03 
^-"^-0.05 


100_6 


0.3210-* 


2+1-7 


13.08 


13 


-1 qn+o-05 


1.20+°-f 


6.4+0-2 


0.0+0.5 


09+0-02 


65+1-0 


2 67+0.14 
^-"'-0.02 




71+0-7 
"■ '1-0.13 


n Sfi+O-54 
O.Sb_(, gg 


10.68 


26 


2 07+°-" 
^■'J' -0.03 


1 07+0.08 
^•^'-0.3 


fi c:+0.25 
'^•"-0.1 


0+0.5 




i.04t[!;i 


9 7O+0.14 

^- '-^-o.os 


100_3 


0-7olo;7 


0.0+1-6 


3.023 



Note. — Parameters values and associated errors, adopting the model DBBEQPAIR for the selected spectra on the HB/NB (first part of the table) and for spectra on the top of the FB 
(second part of the table). The model for the five FB spectra differs from the one adopted for the HB/NB spectra, as it uses a more appropriate disk emission model (DISKPN, in XSPEC) 
instead of the DISKBB component. For the DISKBB component we list the derived apparent inner disk radii in km, while for the spectra fitted with the DISKPN component, the inner disk 
radius has been fixed to 6 Rg and the second column shows the value of the normalizion of the component in units of Mq^ cos(i)/(D2/3'l x 102), where the distance is given in kpc, /3 is 
the color/effective temperature and i is the inclination angle of the disk. The unabsorbed flux of the EQPAIR component has been calculated in the 0.1-200 keV energy range and it is in 
units of 10~* erg cm2 s~i. See Section |4]for an explanation of all the other parameters involved in the fits. All the uncertainties are calculated at 90% confidence level. 



